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The single impurity Anderson model (SIAM) is studied within an enhanced non-crossing approx- 
imation (ENCA). This method is extended to the calculation of susceptibilities and thoroughly 
tested, also in order to prepare applications as a building block for the calculation of susceptibilities 
and phase transitions in correlated lattice systems. A wide range of model parameters, such as 
impurity occupancy, temperature, local Coulomb repulsion and hybridization strength, are studied. 
Results for the spin and charge susceptibilities are presented. By comparing the static quantities 
to exact Bet he ansatz results, it is shown that the description of the magnetic excitations of the 
impurity within the ENCA is excellent, even in situations with large valence fluctuations or vanish- 
ing Coulomb repulsion. The description of the charge susceptibility is quite accurate in situations 
where the singly occupied ionic configuration is the unperturbed ground state; however, it seems to 
overestimate charge fluctuations in the asymmetric model at too low temperatures. The dynamic 
spin excitation spectrum is dominated by the Kondo-screening of the impurity spin through the 
conduction band, i.e. the formation of the local Kondo-singlet. A finite local Coulomb interaction 
U leads to a drastic reduction of the charge response as processes involving the doubly occupied 
impurity state are suppressed. In the asymmetric model, the charge susceptibility is enhanced for 
excitation energies smaller than the Kondo scale Tk due to the influence of valence fluctuations. 

PACS numbers: 75.40.Gb,71.27.+a, 72.10.Fk 



I. INTRODUCTION 

The single impurity Anderson model (SIAM) de- 
scribes an impurity of localized /-states with local 
Coulomb interaction embedded into a metallic host of 
non-interacting c-band electrons.^ In its simplest version 
it discards the possibility of a complex orbital structure 
of the impurity and models the local /-states through a 
two-fold degenerate s-orbital. The Hamiltonian for the 
impurity reads 



U 



(1) 



with (/cr) and = f^fa- the usual creation (annihi- 
lation) and number operators for /-electron with spin a, 
respectively. The local one-particle energy is given by e-^, 
and the local Coulomb interaction is the usual density- 
density interaction proportional to the matrix element U. 
The non-interacting conduction electrons are modeled by 
a single band of Bloch states with crystal momentum k 
characterized by the dispersion relation e 



(2) 



These two parts mix via a hybridization amplitude 14 , 
X>=-^^(li/tcfc<, + /i.c.) . (3) 



The total Hamiltonian is then the sum of these three 
terms 



V 



(4) 



Even though the thermodynamics of the model can 
be solved exactly within the Bethe ansatz methods- 
dynamic quantities can in general not be obtained ex- 
actly and one has to rely on approximations. The 
SIAM has been extensively studied with various meth- 
ods, including the numerical renormalization group 
(NRG),^^^ the (dynamic) density- matrix renormalization 
group ((D-)DMRG)^i^ quantum Monte Carlo (QMC) 
method^i^— and direct perturbation theory with re- 
spect to the hybridization — Especially with the devel- 
opment of the dynamical mean- field theory (DMFT),^^- 
where the solution of an effective SIAM represents the es- 
sential step towards the solution of the correlated lattice 
system, the interest in accurate and manageable impurity 
solvers has increased. 

In this work, we extend the well established enhanced 
non-crossing approximation (ENCA)22i~— to the calcu- 
lation of the static and dynamic susceptibilities of the 
impurity. Like many other approximations formulated 
within the direct perturbation theory with respect to 
the hybridization^i2S — the ENCA is thermodynami- 
cally conserving in the sense of Kadanoff and Baym^^i^ 
It extends the usual non-crossing approximation (NCA) 
to finite values of the Coulomb repulsion U via the in- 
corporation of the lowest order vertex corrections, which 
are necessary to produce the correct Schrieffer- Wolff ex- 
change coupling and the order of magnitude of the low 



2 



energy Kondo scale of the problem. From the NCA 
it is well known that some pathological structure ap- 
pears at the Fermi level below a pathology scale^i^ 
Tpath ^ 10~^ — 10~^Tk- The ENCA removes the cusps in 
spectral functions associated with this pathology^^ and 
only a slight overestimation of the height of the many- 
body resonance at very low temperatures remains. As 
it will be shown in this work, the skeleton diagrams se- 
lected within the ENCA suffer from an imbalance be- 
tween charge and spin excitations and overestimate the 
influence of charge fluctuations. Other than that, it has 
no further limitations. 

Despite the known limitations of the NCA it has been 
widely applied to more complex situations due to the 
forthright possibility of extensions. For the SI AM out of 
equilibrium it is one of the few methods to incorporate 
nonequilibrium dynamics as well as many-body effects.— 
Complex orbital multiplets can be included in a straight- 
forward manner and connections to experimental data 
can be made.^^ However, incorporating finite values of U 
may change the many-body features near the Fermi level 
considerably. ^^^^^ Therfore, a well tested extension of the 
ENCA in order to calculate susceptibilities at finite U 
with the same accuracy as the spectra is desired. In par- 
ticular, the calculation of lattice susceptibilities within 
DMFT— — needs a reliable strategy for an effective im- 
purity, and the treatment of one- and two-particle exci- 
tations of the same footing is of paramount importance. 
Calculations of lattice susceptibilities and phase transi- 
tions of the Hubbard model with the ENCA as the im- 
purity solver are presented elsewhere ^^i^ 

Compared to the "numerically exact" schemes like the 
renormalization group methods (RG), exact diagonal- 
ization (ED) or QMC, the direct perturbation theory 
has its advantages: (i) The approximations are free of 
systematic errors stemming from the discretization of 
the conduction band (RG and ED) or imaginary time 
(QMC). The continuum of band states is kept through- 
out the calculation, and dynamic Green functions are for- 
mulated with continuous energy variables. Thus, there 
are no discretization-artefacts^ and there is no need for 
artificial broadening parameters^i^ or averaging 4^*^ 
(ii) The coupled integral equations for dynamic quan- 
tities, which have to be solved numerically, are formu- 
lated on the real frequency axis, which renders the non- 
trivial numerical analytic continuation of a finite set 
of Fourier coefficients^^ or deconvolution^ unnecessary. 
The continuous-time quantum Monte Carlo approach 
(CTQMC)^^ avoids a systematic Trotter -error (i) but it 
is still plagued with the occurrence of a negative sign- 
problem^^ 

The discretization errors (i) are of special importance 
for self-consistent calculations like the DMFT. There, the 
solution of an impurity model is used to construct a guess 
for the Green function of the lattice which is then used to 
yield a new effective "conduction band" for the impurity 
model. Errors in the treatment of the continuum of band 
states will be propagated by the iterative solution and 



may lead to an incorrect distribution of spectral weight. 

The ENCA can be solved quite effectively on simple 
desktop computers, and is numerically not very demand- 
ing. Spectral functions can be calculated within some 
seconds to minutes while dynamic susceptibilities may 
take up to some hours. Additionally, it contains no free 
parameters and no fine-tuning is necessary. This makes 
it especially interesting for involved lattice calculations. 

II. THEORY 

A. Direct perturbation theory 

In direct perturbation theory with respect to the hy- 
bridization term V the "unperturbed" system is repre- 
sented by the uncoupled {Vk = 0) interacting impurity. 
This is diagonalized by the ionic many-body states \M) 

Hf = ^£^M Xm,m , (5) 

M 

where the operators Xm,m = |^)(^| ^ire projectors on 
the eigenstates \M) and are diagonal versions of the so- 
called ionic transfer (or Hubbard) operators Xm,m' = 
\M){M'\. For a simple s-shell the quantum numbers M 
characterize the empty |0), singly occupied with spin a 
I a) and doubly occupied |2) impurity states with the cor- 
responding unperturbed eigenvalues Eq = 0, E^j = 
and E2 = + /7, respectively. Furthermore, the parti- 
tion function and dynamic Green function are expressed 
in terms of a contour- integration in the complex plane, 

Z = Tre-l'" = ^^^e-l^'Trlz-H]-' (6) 

• [z^ir]-H]~^Bj , 

with if] either a fermionic or bosonic Matsubara fre- 
quency depending on the type of the operators A and B. 
The contour C encircles all singularities of the integrand, 
which are situated on the Im z = and lm{z -\- ir]) = 
axes in a mathematicalley positive sense. Performing the 
trace over the c-electrons first, the reduced /-partition 
function Zj, the /-electron one-body Green function 
= Gr .t and generalized susceptibility Xm w = 

Jcr iJcT ' 

G<r V can be expressed as 

y^M M m' m' 




Po{z) Pa{z + iUn) Ao,a(^, i^n) 
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(10) 



Here I\.m{z^z') = Pm{z)Pm{z') and the ionic propaga- 
tors 



Pm{z) = {M\{[z-H^H,] '),|M) (11) 
^ (12) 



z - Em -^m{z) 



which describe the dynamics of an ionic state |M), are 
introduced. In equation (p!2]) the ionic propagator is ex- 
pressed with the help of the ionic self-energy Em(^)- 
(. . . )c = -^Trc{ e~^^^ . . . ) indicates that the trace is 
to be taken over the c-states only, and Zc represents the 
partition function of the isolated c-system. In equations 
(j9j) and ([TOj), Am,M' and Tm,M' represent vertex func- 
tions to be specified later. These equations are graphi- 
cally represented in Figure [TJ 

After re- writing the Hamiltonian in terms of the ionic 
transfer operators Xm,m', the resolvent operator is ex- 
panded with respect to — i^] ^ = [z — Hf — 



n- 



H, 



Consequently in- 
^.m\M)\c){c\{M\ and using 



V[z-Hf 
serting the identity 1 = XI 

the representation (fTT]) for the ionic propagators, this 
perturbation theory can be formulated with time-ordered 
Goldstone diagrams, representing the dynamics of the 
ionic states \M) and their hybridization processed. 

Approximations are then introduced for the self- 
energies Sm and vertex functions Am,M' and Tm,M'- 
In order to be able to describe non-perturbative many- 
body phenomena like the Kondo effect, certain classes 
of diagrams have to be re-summed up to infinite order, 
resulting in a formulation in terms of skeleton diagrams, 
and corresponding coupled integral equations for the rele- 
vant dynamic functions. Deriving these equations consis- 
tently through functional derivatives from one Luttinger- 
Ward functional <l> renders these approximations thermo- 
dynamically consistent. 



B. ENCA for generalized dynamic susceptibilities 



Within the ENCA, the ionic self-energies and propa- 
gators have to be determined from the well known set of 



coupled integral equations ^i^^i^^i^ 

^o(^) = ^ fdxA-^{x)P^{z^x)AoA^^^) (13) 
Ecr(^) = J dx A~{x) Po{z — x) Ao,cr(^ — ^, x) 

^2{z) = Y^ dxA-{x)Pa{z-x)A2,a{z,x) , 
cr 

where the vertex functions are given by 

Ao,a(^, ^0 = 1 + Jdx A+(x) P^{z + x) P2{z + + x) 

(14) 

A2,a(^,^0 = l+ dxA-{x)P^{z-x)Po{z-z-x) , 



and A^{x) = A(x)/(±x). The hybridization function 
A{x) is constructed from the c-band electrons 



(15) 



and f{x) = l/(e^^ + 1) is the Fermi function with /3 = 
1/T the inverse temperature (note h = c = ks = l)- 

For the generalized susceptibilities Xm,m'{z, z')^ an ad- 
ditional set of integral equations has to be set up which 
is shown graphically in Figure [2l In principle, there 
are 16 such functions, but for each quantum number 
M only the four functions Xm,q{z^ z')^ XM,a{z^ z')^ and 
Xm,2{z^z') are coupled. Due to the conserving nature 
of the ENCA, these equations are closely related to the 
ENCA expressions for the ionic self-energies; only some 
additional bosonic lines entering and leaving the site have 
to be introduced. 

Whereas the lowest order vertex corrections of equa- 
tion (fT4|) together with the ionic propagators (pTj) and 
self-energies (p!2)) . are sufficient to furnish a conserving 
approximation, the vertices T m,M'{z^ z') for the suscep- 
tibilities ([To]) have to be iterated up to infinite order! The 
transcription of these graphs into formulas is straight for- 
ward but lengthy and will be omitted here. 

For very large Coulomb interactions the spectral 
weight in the doubly-occupied propagator is completely 
moved to energies of the order of U . This leads to a neg- 
ligible infiuence in the region of accessible energies of the 
order of the bandwidth due to the large energy denom- 
inators, and all terms involving P2 effectively drop out 
of the equations. Therefore, for infinitely large Coulomb 
repulsion ^ 00, the ENCA reduces to the usual NCA 
and all equations presented above approach the ones al- 
ready known from the literature ^^i^ In this sense, the 
ENCA is still referred to as non-crossing even though 
crossing diagrams are included. For approximations in- 
volving "real" crossing diagrams see Refs. l23ii29ll3Qll53l . 
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FIG. 1: Diagrammatic equations for the ionic self-energies (upper), /-electron Green function (lower left) and dynamic suscep- 
tibility (lower right, Xo,a{ii^) as an example). The dashed, wiggly and double- wiggly lines represent the empty, singly occupied 
{a represents the opposite of the spin a) and doubly occupied ionic propagators, respectively. The triangles and boxes represent 
the appropriate vertex functions. The closed full lines in the upper graphs represent uncorrelated propagations in the c-band. 



The functions Xm,M'{z^ z') fulfill the symmetry rela- 
tions 



and 



Xm,m'[z ,z ) = Xm,m'[z,z) 



(16) 



(17) 



which are revealed by inspection of the perturbation ex- 
pansion. The susceptibilities also obey the sum rules 



M' 



(18) 
(19) 



M 



XM' .M^'^^n) = ^^^{Xmm) 5^Vr^,Q = ^^Vr^,Q •> (20) 



M 



which arise from the completeness relation of the lo- 
cal ionic states, 1/ = XIm These sum rules 
transform into equivalent statements for the functions 
Xm',m{z, z'), 



Pm{z)-Pm{z') 



(21) 



M' 



It can be analytically checked that the ENCA respects 
these sum rules. The form of equation ([2T]) explicitly 
reveals the conserving nature of this approximation, since 
the sum of sets of two-particle correlation functions is 
determined by the corresponding one-particle correlation 
function, i.e. the ionic propagator. Insofar the equations 
(|2T]) resemble generalized Ward-Identities. 

In order to obtain the physical susceptibility as a func- 
tion of a real frequency, the contour integration of equa- 
tion (p!Q]) has to be performed, and then the external 



bosonic Matsubara frequency can be analytically contin- 
ued to the real axis, ivn v ^ iS = (^ > infinitesi- 
mal). The result of this procedure reads 



p du 



lM,M'(c^,a; + i/) (22) 



where the symmetry relations (jT6|) were used and the 
defect quantities 



YM,M'i^^y) 



-(3x 



XM,M'i^^^y^) - XM,M'ix ,^+) 

(23) 



were defined. Due to the appearance of the exponential 
factor, a direct numerical evaluation of the defect quan- 
tities given the Xm,M' possible, and additional sets 
of integral equations have to be solved for the Ymm' (cf. 
Ref.lEiandlHi). 

In the following, only spin-symmetric situations are 
considered, and the propagators for opposite spins are 
identified, i.e. = = P^, Xa.a = Xa,a Xt,t' Xa,a 

X(T,a = Xt,i 



1. Magnetic susceptibility 

The relevant quantity for the magnetic susceptibility is 
the auto-correlation function of the z-component of the 
spin operator ^ n^ — fii = X^^^—X^^i. This translates 
into the linear combination 



Xmag{z, z') = Xt,t(^. ^0 " Xt,i(^. ^0 



(24) 



FIG. 2: System of integral equations for xm,o{z, z^), XM,a{z, z')^ and xm,2(^, z'^ within the ENCA. Zig-zag lines in the vertex 
part represent the ionic state M, which can either be empty, singly occupied or doubly occupied. The arrows entering and 
leaving each diagram carry the external bosonic frequency iv. 



which needs to be determined. Setting up the equations The derivation of the corresponding defect equation for 
for this linear combination using the general equations Ymag along the lines of the definition ([23|) is straight for- 
depicted in Figure [2] leads to the cancellation of all spin- ward, but will be omitted here for brevity, 
symmetric terms. The function Xma^(^,^0 decouples 
from the Xo,m(^,^0 X2,M^^^^'\ and only one inte- 
gral equation results, 




(25) 



For the charge susceptibility the relevant quantity is 
the auto-correlation function of the charge operator h = 
+ — ^t,t + + 2X2,2 leading to the dynamic 



2. Charge susceptibility 



+ Po{z-y) P2{z^x) Xmag{z^x-y,z'^x-y) 



6 



function 



Xchargeiz, z') = Xt,t(^' ^0 + Xt,i(^. ^0 (26) 

= X2,2(^, ^0 - X0,2(^, Z') + XO,o(^, ^0 - X2,o(^, Z') 
+ ^O,t,2(^,^0 . 

In the second equality sum rules like ([2T]) were used. The 
function £'0,^,2(^5 ^0 is not relevant for dynamic suscepti- 
bilities, since it only contributes in the case of a vanishing 
external frequency iun = 0. 

The hole susceptibility, i.e. the auto-correlation func- 
tion of the hole operator h = 2 — n = 2Xo,o + -^t,t + ^i,i' 
is given by the same expression, only with a different con- 
tribution from the sum rules. Both of these contributions 
do not change the dynamic susceptibility since they pro- 
duce only the necessary constant shift between the static 
{ivji = 0) charge and hole susceptibility after the contour 
integration, 

Xholeii^n) = Xcharge{i^n) + 4(1 - (n))(5i^^,o • (27) 

The dynamical quantities of interest are best obtained 
by setting up the integral equations for the linear com- 
binations co,o = Xo,o - X2,o, co,t = Xo,t - X2,t and 
Co, 2 = Xo,2 - X2,2, which read 



-f Ao,t(^,:^^) - 1 co^^{z^x,z^x) (28) 
^2 J dxdyA-^{x)A-^{y) P^{z'^x) P^{z^y)x 

X co,2{z^x^y,z'^x^y) > 



co,t(^,^0=nt(^,^0| (29) 

JdxA~{x) Ao^^{z^ — x^x) 

+ Ao^^{z — x^ x) — 1 co^o{z — x^ z' — x) 

+ JdxA'^{x) [A2^^{z ^x,x) 

+ A2,t(^+a:, x) - l] €0,2(^+2;, z'-^x) 
-h jdxdyA'^{x)A~{y) P2{z' ^x) PQ{z-y) 

Po{z'-y) P2{z^x) co^^{z^x-y,z'^x-y) > 



co,2(^,^0 = n2(2;,2;0| -1 + 2 J dxA' (x) ^A2,t(^',^) 

+ A2,t(^,a:) - 1 co^^{z-x,z-x) (30) 
^2 J dxdyA-{x)A-{y) P^{z-x) P^{z-y) x 

X co^o{z-x-y,z-x-y) } . 



All three linear combinations are now coupled which 
makes the numerical solution of the full system in- 
evitable. 

The corresponding defect equations are again deter- 
mined in a straight-forward way, but are omitted here 
due to their length. 



C. Numerical evaluation 

The ionic propagators and defect propagators are 
strongly peaked at the ionic threshold, and they even 
develop a singularity at zero temperature j^i^i^ The po- 
sition of this ionic threshold is not know exactly a priori, 
so we use self-generating adaptive frequency meshes to 
represent all functions numerically. 

After obtaining the ionic propagators from the set 
of integral equations Eq. ([nj), Eqs. (j25j) and (j28]) - 
(f3Q|) can be solved to yield the two-particle quantities 



Xmac^(^^7 ^ i ^^) and cm,M'(^^, ^ i ^^)- With these at 
hand, the corresponding defect equations can be solved 
and physical susceptibilities along the lines of Eq. (|22]) 
can be extracted. 

All integral equations are solved via Krylov subspace 
methods^ where — starting from a sophisticated guess 
for the unknown functions — the equations are iter- 
ated until convergence is reached. In order to ac- 
celerate convergence and stabilize the iteration proce- 
dure Pulay-mixing^ between different iterations is im- 
plemented (see, for example, Ref . IHol or [6Q|) . 

Unfortunately, the system for the defect quantities be- 
comes singular for vanishing external frequency u ^ 0. 
This is seen at the sum rule (|2T]) . which translates into 
J2m' ^M',m(^,^ + z^) = -27ri^M{(^)/{^ + iS) implying 
that the Ym,m' become of the order of ^ for small 
while all terms in the integral equations stay at the order 
one. This becomes of some importance when extracting 
static susceptibilities from calculations of the dynamic 
susceptibilities, where a very small but finite frequency 
is used. The resulting convergence problems are proba- 
bly connected to the ones already mentioned in Ref. |6ll. 
More details on the numerics can be found in Ref. |4l|. 
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FIG. 3: Rescaled /-electron spectra (upper graph) and total 
self-energy (lower graph) as a function of energy for various 
values of the Coulomb repulsion U. The temperature was 
fixed at half the Kondo temperature T ^ Tk/^ of each case, 
the ionic level was chosen e-^ — —1 and the Anderson width 
was Aa = 0.2. The conduction band DOS was that of a 
three-dimensional tight binding lattice with half bandwidth 
W = 3 ("c-band" in the upper graph). The insets show the 
low energy interval around the Fermi level, where the energy 
is measured in units of the corresponding Kondo temperature 
and the horizontal lines indicate the T = limiting values. 



III. RESULTS 

As it is well known, the SIAM exhibits the Kondo ef- 
fect for sufficiently large Coulomb repulsion U and the 
single particle level below the Fermi energy, < 0. A 
very prominent manifestation of that effect is the emer- 
gence of the Abrikosov-Suhl resonance (ASR) in the one- 
particle density of states (DOS) at the Fermi level for 
temperatures of the order of the Kondo temperature Tk- 
Within the ENCA an order of magnitude estimation for 
Tk is given hy^ 



(31) 



where J is the effective antiferromagnetic Schrieffer- Wolff 
exchange coupling, the Anderson width = 7rpc{0)V'^ 



and W the half bandwidth of the conduction electron 
DOS pc{^)- The hybridization matrix element was as- 
sumed to be local, i.e. momentum independent Vk^V. 

As already mentioned, the pathology of the NCA mani- 
fests itself in the overestimation of the height of the ASR 
and a violation of Fermi liquid properties for too low 
temperatures as well as in in situations with large va- 
lence fluctuations j^i^i^i^ This pathological behavior is 
strongest for the case of a spin-only degeneracy {N = 2), 
which is considered in this work. 

Figure [3] shows the one-particle /-electron DOS 
pf{uj) = ^lmF(j{uj — iS) and the imaginary part of the 
total self-energy ImT^y* {uj — iS) = —lm[l/F(j{oJ — iS)] cal- 
culated within the ENCA and the NCA {U = oo). As 
can be seen, the ENCA does not overestimate the height 
of the Kondo peak or violate the Fermi liquid property of 
the total self-energy ImSj'^(0 — iS) > A a down to tem- 
peratures of half the Kondo temperature. In contrast, 
the U = oo NCA curves in the graphs do violate these 
limits for the same parameter values. 

The fact that the ENCA performs better than the NCA 
in comparable situations is a consequence of a better bal- 
ance between different kinds of perturbational processes. 

The importance of such a balance has been pointed out 
repeatedly.23^62,63 

Even though the performance of the ENCA is consider- 
ably enhanced over the NCA, it eventually overestimates 
the height of the ASR for even lower temperatures and 
still misses to produce the correct T = Fermi liquid 
relations. Further improvements can be archived via the 
incorporation of higher-order diagrams 



A. Benchmarking the ENCA 

1. Static susceptibilities in the symmetric case 

In order to obtain a better understanding of possible 
shortcomings of the ENCA it is worthwhile to consider 
charge and spin excitations separately and benchmark 
them against some exactly known results. The thermo- 
dynamics of the SIAM can be obtained exactly from the 
Bethe ansatz method^^i^i^— At zero temperature and 
for the symmetric case {e^ = — with a flat conduction 
band of infinite bandwidth {W oo) the static suscep- 
tibilities can even be obtained in closed form^^iS 



exact frri r.\ 
imag U - - ^ I 1 



1 



3 16x } 



2U g I 

dx — 



with 



exp < 



Till TtA/ 



8A^ 



2U 



(32) 



(33) 
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FIG. 4: Static magnetic (upper) and charge (lower) suscepti- 
bility for the symmetric case e-^ = — [7/2 as a function of U for 
finite temperatures Tk and Tk/2 in comparison to the exact 
T = results of equations (|32)) and (|34|) . The calculations 
were done for a constant c-electron DOS with half bandwidth 
W = 10 and an Anderson width Aa = 0.3. 

and 

(34) 

Apart from the small coupling correction ^ ^ in the 
exponent, Tl exactly coincides with the Kondo temper- 
ature of equation (|3T]) . 

In the graphs of Figure 2] the exact zero temperature 
magnetic (upper) and charge (lower) susceptibilities of 
equations (|32]) and (solid red lines) are compared to 
the ENCA susceptibilities (blue dots) for e-^ = — [7/2 as 
a function of U . For the ENCA curves two characteristic 
temperatures Tk and are chosen. 

The characteristic exponential U dependence of the 
magnetic susceptibility is essentially the same as for the 
exact Bethe ansatz result, but the absolute height is 
somewhat different. At T = Tk the magnetic suscep- 
tibility is roughly 70% of its zero temperature saturation 
value. However, at T = Tx/2 the susceptibility has al- 
most saturated and the agreement with the T = 0- value 
is quite good. 

The deviations for /7 > 4 are due to the method used 
to extract the static susceptibilities: The static limit is 
obtained by evaluating the dynamic susceptibility at a 
small but finite external frequency, in our case v = 10~^. 
Since the ENCA represents a conserving approximation, 
the results obtained with this method agree with the 
static susceptibility obtained from a derivative of a ther- 
modynamic potential, or solving separate equations as in 
Otsuki et al.— . This is valid as long as the minimal fre- 
quency is negligible compared to the lowest energy scale 
in the problem. 

However, the Kondo temperature for = 4 is only 



FIG. 5: Magnetic and charge susceptibility for the symmetric 
impurity {e^ — —U/2) as a function of the Coulomb interac- 
tion [/ at a fixed temperature T = 0.05 and for an Anderson 
width Aa = 0.2. Notice, the x-axis is scaled as U/{1 + U), 
so that U — oo corresponds to U/(l -\- U) = 1. The arrow 
on the right of the figure corresponds to the value of a Curie 
susceptibility Xmag / {qij^bY — 1/(4^) = 5 of a free spin at 
temperature T = 0.05. 

of the order of 10~^ and decreases for larger U. The 
susceptibility calculated at = 10~^ then does no longer 
represent the static limit anymore, and the decrease seen 
in Figure m is produced, which is therefore not indicating 
a shortcoming of the ENCA method. The deviation can 
be cured by choosing a smaller value for the external 
frequency. 

The charge susceptibility (lower graph in Figure 2]) 
shows no significant temperature dependence for T = Tk 
and Tk/"^ and lies right on top of the exact T = result. 
The deviations for /7 > 4 are explained in the same way 
as for the magnetic susceptibility described above. 

These reasons for the deviation at larger U can be fur- 
ther confirmed by calculating the static susceptibilities 
for a fixed finite temperature T = 0.05 as a function of 
the Coulomb interaction, which are shown in Figure [5l 
The figure compares the static magnetic (red dots) and 
charge (blue dots) susceptibilities for the symmetric sit- 
uation with the exact T = results. 

The charge susceptibility decreases monotonically with 
increasing Coulomb repulsion as expected. It follows 
the exact zero temperature susceptibility very accurately, 
which again confirms its weak temperature dependence 
in symmetric situations at low temperatures. 

The magnetic susceptibility (red dots) does agree with 
the exact T = solution for low Coulomb repulsions 
but deviates for U/{1 + /7) > 1/3 (U > 1/2). This is 
a finite temperature effect, since for values of U > 1/2 
the Kondo temperature is smaller than the chosen finite 
temperature of T = 0.05. Consequentlya, for U > 1/2 
we are not in the low temperature regime, and the sus- 
ceptibility is not well described by its T = value. The 
magnetic susceptibility does not grow exponentially with 
U as for T = 0, but instead saturates for U ^ oo at 
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FIG. 6: Temperature dependent squared effective magnetic 
moment within the ENCA (colored dots) for a fixed Coulomb 
repulsion U — AV'^ and three different /-level positions 
= 0, — V^^, —2V'^. The solid grey curves are the exact Bethe 
ansatz results for the same parameter values (after Okiji and 
Kawakami^^). The calculations were done for a constant c- 
electron DOS with half bandwidth Vl^ = 10 and A a = 0.167. 



FIG. 7: Static magnetic susceptibility for a fixed T = 0.05 
and A A = 0.2 as functions of the ionic level position e-^ rela- 
tive to the half- filling value —U/2 for various values of U. The 
conduction band was chosen to be constant with a half band- 
width of = 10. Curves without dots ("free") are calculated 
without two-particle interactions, i.e. with the particle-hole 
propagator of equation (|36|) . 



the asymptotic value of the Curie susceptibility of a free 
spin, Xmag = 1/ (4T) = 5, which is indicated by an arrow 
on the right border of the graph. 

2. Static susceptibilities in the asymmetric case 

As a first test for the ENCA in the asymmetric sit- 
uation Figure [6] compares the square of the effective 
screened local magnetic moment of the impurity, 

to the exact Bethe ansatz result for three different ionic 
level positions as a function of temperature. The solid 
grey lines are the exact Bethe ansatz solution, which is 
taken from Ref. [69|, while the colored points are ENCA 
calculations for the same parameter values. The half 
bandwidth was taken to be W = 10, which should be 
large enough to be comparable to the Bethe ansatz solu- 
tion where = oo. 

The ENCA slightly overestimates the squared effec- 
tive moment but all characteristic features are essentially 
the same as for the Bethe ansatz. Especially the shape 
and the relative height of the curves is in remarkable 
agreement: All three ENCA curves can be brought to lie 
right on top of the exact Bethe ansatz results when they 
are rescaled with one single factor. This indicates that 
the ENCA produces a slightly modified Kondo scale, but 
otherwise describes the static magnetic properties almost 
exactly. This is especially remarkable for the intermedi- 
ate valence situation with e-^ = 0, where the empty and 
singly occupied ionic configurations are almost degener- 
ate. In such situations stronger pathologies occur in the 



one-particle DOS and the NCA-type of approximations 
would be expected to yield results of lower quality. How- 
ever, as it can be seen, the magnetic excitations are still 
described very accurately. 

Calculations of the effective squared local moments for 
a wide range of Coulomb interactions and hybridization 
strength (not shown) resemble the exact solutions known 
from the literature^i^i^ and the characteristics of the 
different asymptotic regimes are well reproduced by the 
ENCA. 

In Figure [71 the magnetic susceptibility is systemati- 
cally examined as a function of the ionic level position 
6^ and for fixed T and U. Also shown are the suscep- 
tibilities without explicit two-particle interactions (lines 
without dots labeled as "free"), i.e. the local particle- hole 
propagator P'^(O) calculated via 

/oo 
dco f{oj)pf{co) (36) 
-OO 

X [F^{uJ + + i(5) + F^{u - V - id)] , 

where represents the one-particle /-Green function 
(cf. equation (|9])), pfioo) the corresponding spectrum and 
fiuo) the Fermi function. 

All curves are symmetric around \U = which 
just reflects the particle- hole symmetry of the model. 
The particle-hole propagators shows the expected max- 
ima approximately situated at the positions of the Hub- 
bard peaks in the one-particle DOS ^ \U ~ 

For very small U = 0.01, the susceptibility calculated 
with the ENCA is indistinguishable from the particle- 
hole propagator, as a consequence of the near lack of 
two-particle correlations. 
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FIG. 8: lim.^o o ^'^""f^^^ ^"^^^"^^"'^^ evaluated as a func- 
tion of temperature and various Coulomb repulsions within 
the ENCA. For comparison the U = oo-NCA curve is shown 
as well. A 3d-SC DOS was used for the conduction electrons 
and Aa = 0-3. 

For larger Coulomb interactions, the ENCA suscepti- 
bility shows only one broad maximum around -\-^U = 
(half filling), which grows in height and width with in- 
creasing U. The enhancement of the susceptibility is due 
to the increasing local magnetic moment with larger U. 
The plateau which develops around zero, is due to the 
stability of the local moment as long as the singly occu- 
pied ionic configuration is stable, i.e. the lower Hubbard 
peak being below and the upper above the Fermi level, 
and the temperature is not too low compared to Tk- But 
as soon as one of the Hubbard peaks extends over the 
Fermi level, i.e. + |Aa <0 ot - |Aa > 0, the 

moment is destabilized. For both Hubbard peaks below 
(above) the Fermi level, the impurity is predominantly 
doubly occupied (empty) and the magnetic susceptibil- 
ity drops drastically. The curve then rapidly approaches 
the particle- hole propagator, indicating that explicit two- 
particle interactions are unimportant. 

The reproduction of the correct results for the effec- 
tively non-interacting limit at U = 0.01 as well as the 
empty- or fully occupied regimes is quite remarkable. 

As it was already mentioned earlier, the NCA does 
violate Fermi liquid properties for very low temperatures. 
Another indication, in addition to the imbalance of the 
imaginary part of the total self energy at the Fermi level 
(-ImE(0 + i(5) < A^), stems from the zero frequency limit 
of the imaginary part of the dynamic susceptibility. For 
the Fermi liquid at T = it has to obey the so called 
Korringa-Shiba relation^ 

hm is^B? lmxma,{^) 

The function on the left hand side of this relation is shown 
in Figure [8] as a function of temperature for various values 
of U . The explicit form ([37]) holds for a flat infinitely- 
wide conduction band DOS. For a different DOS the nu- 
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FIG. 9: Temperature dependent static charge susceptibil- 
ity XchargeiS^) for various Coulomb repulsions U (Aa = 0-2, 
= —1). The arrow at the left border shows the exact T = 
limit for the symmetric case U — 2. The high temperature 
asymptotic form of a non-interacting impurity (x = 1/(8T)) 
is shown as well (thick line). The colored dots, which are 
not connected by a line, denote the asymptotic susceptibili- 
ties calculated without any two-particle interactions, i.e. from 
P-^(zy = 0), see equation (|36)) . The conduction band was cho- 
sen to be that of a 3d-SC lattice with a half bandwidth of 
W = 3. 



merical pre-factors might change slightly, but the left- 
hand side is still expected to be of the order of unity due 
to universality of the SI AM at low energies. 

For the NCA, the quantity \imj^^olmXmag{^)/^ is 
known to diverge at T = 0^ which is reproduced by 
the U = oo curve in the figure. However, the ENCA 
(finite U values) performs considerably better than the 
NCA. The curves still slightly increase for temperatures 
T ^ Tk/2^ but they eventually saturate at a finite value 
and do not diverge. This represents a considerable im- 
provement of the qualitative behavior of the ENCA over 
the NCA. 

The temperature dependent static charge susceptibil- 
ity is shown in the Figure [9] for various values of U as 
a function of temperature. For high temperatures and 
U < A the static susceptibility behaves effectively non- 
interacting with Xcharge = 1/ (8^) ^S CXpCCtcd. 

For U = 4 the susceptibility still has the characteristic 
1/T dependence for high T, but with a prefactor more 
closely to 1/16. This can be understood since in this sit- 
uation the upper Hubbard peak (incorporating roughly 
half the spectral weight) is energetically just above the 
upper band edge of the c-band and therfore the accessi- 
ble spectral weight for two-particle excitations is approx- 
imately halved. 

The rapid drop of the susceptibility for U = 4 at tem- 
peratures T > W is attributed to inaccuracies in the 
numerics for solving the integral equations. However, in 
this effectively non-interacting regime the susceptibility 
can be calculated without explicit two-particle interac- 
tions via the particle- hole propagator of equation ([36|) . 
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The results thus obtained are shown in the graph as col- 
ored dots without joining lines and are seen to be nicely 
proportional to 1/T for high temperatures. 

Therefore, the ENCA nicely reproduces the high tem- 
perature asymptotics of the SI AM. 

At temperatures around T/W ~ | all susceptibil- 
ities show a pronounced maximum, which stems from 
thermally excited charge fluctuations between the empty 
and singly occupied ionic levels with excitation energy 
^ \e^/W = 1/3. For U = 2, fluctuations between the 
singly and doubly occupied states have the same excita- 
tion energy |e-^|/W^ = (e-^ + U)/W = 1/3 and therefore 
contribute equally. For [7 = 1.7 the energy of fluctuations 
involving the doubly occupied state is somewhat smaller 
M ~ 0.23) and the peak is therefore broadened to 
lower energies. For [7 = 4 the doubly occupied state is 
inaccessible for thermal fluctuations; so only the empty 
and singly occupied levels contribute leading to a reduc- 
tion of the susceptibility maximum by approximately a 
factor of two. 

At lower temperatures {T/W < 0.05 in the figures) the 
charge susceptibilities exhibit a slow increase followed by 
a saturation at the zero temperature values. The in- 
crease in the charge susceptibility occurs in a tempera- 
ture range, where the Kondo singlet and the local Fermi 
liquid formation take place, which manifests itself in the 
growing many-body resonance at the Fermi level in one- 
particle DOS pfioo). 

Even though a direct interpretation of the increase in 
terms of a Fermi liquid picture (where the charge suscep- 
tibility is proportional to the DOS at the Fermi level) is 
not applicable since the Fermi liquid is formed only at 
very low temperatures, it still provides an intuitive way 
of understanding: The increasing spectral weight at the 
Fermi level leads to an enlarged phase space volume for 
two-particle excitations and the charge susceptibility is 
at least roughly proportional to the DOS at the Fermi 
level. This is supported by the fact, that Xcharge{^) 
increases logarithmically with decreasing temperature, 
which is also the case for p/(0). But how strong the 
increase actually is and how it is influenced by the value 
of the Coulomb repulsion U cannot be deduced from the 
simplified Fermi liquid analogy. This rather depends on 
the details of the two-particle correlations. 

In the symmetric situation (U = —2e^ = 2) the charge 
susceptibility increases only moderately and approaches 
the exact T = limiting value known from the Bethe 
ansatz, which is indicated by the arrow at the left border 
of the Figure [9l For the asymmetric cases, the increase 
is considerably more pronounced. Especially, the dras- 
tic low temperature increase for /7 = 4 is rather unex- 
pected. The absolute value of the susceptibility is even 
larger than for the smaller values of [/, which is counter- 
intuitive since charge fluctuations should be suppressed 
for larger U . However, the tendency that for a given 
level position the charge susceptibility in the asym- 
metric situation can increase with growing U is known 
from perturbation theory^ as a characteristic feature of 
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FIG. 10: Static charge susceptibility for a fixed T = 0.05 and 
A A = 0.2 as functions of the ionic level position relative 
to the half- filling value —U 12 for various values of U . The 
conduction band was chosen to be constant with a half band- 
width oiW — 10. Curves without dots ("free") are calculated 
without two-particle interactions, i.e. with the particle-hole 
propagator of equation (|36p . 



valence fluctuation physics. Valence fluctuations being 
at the origin of this enhanced low temperature increase 
of the charge susceptibility are in agreement with the 
observation already made above, that for /7 = 4 the dou- 
bly occupied ionic orbital is outside the conduction band 
and the system is therefore from the outset closer to the 
intermediate valence fixed point. 

Reference calculations with the NRG (not shown) in- 
deed display the characteristic features of the charge sus- 
ceptibility as shown in Figure [9l A maximum for tem- 
peratures of the order of the ionic level positions |e-^| 
and -\-U and an increase towards lower temperatures. 
In situations close to the valence fluctuation regime this 
increase leads to an enhancement of the charge suscepti- 
bility by a factor of about 10. However, the parameter 
values chosen for the U = 4 ENCA-curve are not very 
close to the valence fluctuation fixed point. This is also 
reflected in the magnetic susceptibility for U = 4, which 
does not show any signatures of the valence fluctuation 
regime, but rather exhibits behavior characteristic for the 
transition from a local moment to the strong coupling 
fixed point (not shown). NRG calculations with param- 
eter values similar to the ones chosen in this study, did 
show a low temperature increase, but not as strong as 
observed with the ENCA. 

Altogether it can be concluded, that the ENCA does 
describe the charge fluctuations qualitatively right, but 
overestimates the influence of intermediate valence phe- 
nomena at very low temperatures in the asymmetric case. 

To make the range of applicability of the ENCA more 
clear, it is instructive to consider the charge susceptibility 
for fixed values of U and T, varying the ionic level posi- 
tions e-^, which is shown in Figure [TOl The particle- hole 
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propagators already displayed in Figure [71 are included 
as well (lines without dots). The ENCA charge suscepti- 
bilities are always minimal for half filling (e-^ + /7/2 = 0) 
and increases away from the symmetric case. The abso- 
lute value of the charge susceptibility in the symmetric 
situation is drastically reduced compared to the corre- 
sponding particle-hole propagator for large values U = 2 
and /7 = 4, which indicates, that the two-particle corre- 
lations strongly suppress charge fluctuations. In that sit- 
uation, the susceptibility cannot accurately be described 
by the one-particle DOS alone and independent though 
strongly renormalized quasiparticles. 

On the logarithmic scale, the increase with growing 
distance from zero can nicely be fitted with a parabola 
centered at zero, which corresponds to an exponential 
increase of the susceptibility, Xcharge ^ e"^^ ^ (^q^ > 

0). This shows the strong influence of the asymmetry 
and the contribution of valence fluctuations to the charge 
fluctuations. 

However, the ENCA clearly fails for large asymmetries 
as the susceptibility saturates for -\- ^U\ > ^U^^. In 
contrast, Xcharge{^) should decrease again (cf. Kawakami 
and Okiji— ) and approach the particle- hole propagator, 
due to the effective non-interacting nature. This is most 
drastic for the almost non-interacting case with U = 0.01, 
where, apart from reproducing the value at half filling 
quite accurately, the curve goes the opposite direction as 
expected. The values at which the downturn in the sus- 
ceptibility should occur correspond to situations, where 
both Hubbard peaks in the one-particle spectrum (very 
roughly at and -\-U) are either below or above the 
Fermi level, corresponding to the empty- and fully occu- 
pied impurity regimes. 

The ENCA is designed to describe spin flip scatter- 
ing and the magnetic exchange coupling correctly but 
it does not fully capture the physics of charge fluctua- 
tions outside the Kondo regime. In situations, where the 
unperturbed ground state is either the empty or dou- 
bly occupied ionic state, crossing diagrams neglected in 
the ENCA are vital to describe charge fluctuations accu- 
rately. On the other hand, magnetic fluctuations are still 
described very accurately in these situations (see Figures 
ElandEj). 



B. Dynamic susceptibilities 

1. Magnetic susceptibility 

The imaginary part of the dynamic magnetic suscep- 
tibility is shown in Figure [TT] for two different values of 
U and two characteristic temperatures. The spectrum of 
the susceptibility shows a pronounced maximum, which 
is shifted to lower frequencies and increases considerably 
in height as the temperature is lowered. For tempera- 
tures below the Kondo temperature, the position of the 
maximum remains fixed at a value of the order of the 
Kondo temperature. 
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FIG. 11: The imaginary part of the dynamic magnetic sus- 



ceptibility for 



-1 two values of [/ (t/ = 2, 4) and two 



characteristic temperatures (T = 10Tk,Tk/2) in a double- 
logarithmic plot. The corresponding Kondo temperatures are 
indicated as arrows on the frequency axis. All curves are cal- 
culated for A A = 0.3 and a 3d-SC band DOS (H^ = 3). 



Also shown in the figure are fits with a Lorentzian form 



Xmag' (^) 



Xo 



1 - iv/V 



(38) 



which describe the low frequency susceptibilities very 
well. The form (|38]) corresponds to an exponential spin 
relaxation with relaxation time l/F. The line- width F is 
directly proportional to the NMR impurity nuclear spin- 
lattice relaxation rate F ~ TiJ^ 

The relaxation rates F extracted from susceptibilities 
for various parameters follows a VT-lsiw^ for high tem- 
peratures and saturates at a value of the order of Tk at 
temperatures below Tk (not shown) , in accord with what 
was already found earlier ^^i^i^i 

The physical picture behind these findings is quite 
clear: Upon lowering the temperature, the local moment 
of the impurity becomes increasingly coupled to the sur- 
rounding spin of the band electrons resulting in an en- 
hanced response. At temperatures of the order of or lower 
than the Kondo temperature, the local Fermi liquid state 
is approached in which the local spin is screened and a 
local Kondo singlet is formed with a "binding energy" of 
about Tk' Therefore the maximum in the spin excita- 
tions spectrum, as well as the NMR relaxation rate, both 
are pinned at an energy of the order of Tk- 

Jarrell et al.— also found that the function 



2Xmag(0) 1^ 



(39) 



shows universality and depends only on T/Tk- Figure [12] 
shows this function normalized to its zero frequency value 
for various parameter sets. All graphs can be collapsed 
onto one single curve showing the universal shape of the 
function / for low energies. 
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FIG. 12: Scaling function f{v/T^)/f{id) from equation ^ 
for the following parameter sets: (e-^ = — 0.5;[/ = l; T /Tk — 
10,1), {e^ = -1;U = 2;T/Tk = 1), {e^ = -1.5;U = 3;T/Tk = 
1), {e^ = -2-U = 4;T/Tk = 1), {e^ = -1-U = 3-T/Tk = I), 
(e^ = = 4; T/Tk = 1, 1/2, 1/3, 1/5, 1/7, 1/10) and (e^ = 

-1;[/ = 8;T/Tk = 1). All curves are calculated for Aa = 0.3 
and a 3d-SC band DOS. 



In order to achieve scaling a guess for the actual Kondo 
temperature has to be used. In contrast, the T^- value 
calculated by equation (|3T]) and used in this work does 
not represent the exact physical low energy scale T^, 
but only provides an order of magnitude estimate. In 
the universal regime with a flat c-band DOS this should 
not make any difference, but since we are using the 3d-SC 
DOS, no n- universal corrections enter for different and 
U. The value of the "real" Kondo temperature could have 
been extracted from fits of the calculated susceptibilities 
to the universal curve of the susceptibility as described 
by Jarrell et al. . 

The rapid decrease of the curves in the figure for fre- 
quencies of i^/T^ ^100 also stems from the finite band- 
width of the 3d-SC conduction band used for these cal- 
culations. 

The above findings clearly confirm, that the dynamics 
of the impurity spin is solely determined by the antifer- 
romagnetic exchange between the impurity- and conduc- 
tion electron spins. Even in the asymmetric situation, 
the only relevant energy scale for magnetic fluctuations 
of the impurity is the Kondo temperature Tk at low tem- 
peratures. 



2. Charge susceptibility 

The imaginary part of the dynamic charge susceptibil- 
ities for two different Coulomb repulsions and character- 
istic temperatures, calculated with a 3d-SC band DOS, 
are shown in Figure [131 In the spectra the characteristic 
features stemming from excitations involving the Hub- 
bard peaks at energies around |e-^| and -\-U are clearly 
visible. 

In the symmetric case {U = 2), the height of the peak 
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FIG. 13: Imaginary part of the dynamic charge susceptibility 
for the symmetric (e 



-1, U = 2, upper) and asymmetric 
[e' = — 1, [/ = 4, lower) case for two characteristic tempera- 
tures T = IOTk and T = Tk/2. The curves for the particle- 
hole propagators ("free") are shown as well. The insets show 
the corresponding low energy region on a double-logarithmic 
scale. The Kondo temperature for each case is indicated by 
the red arrow on the frequency axis in the low energy insets. 
The calculations are done for a 3d-SC band and A a = 0.3. 



at ^ |e-^| is about twice of the one in the asymmetric 
case {U = 4). This is due to the doubled phase-space 
volume for the symmetric situation with excitation ener- 
gies matching, -\- U = |e-^|, while for /7 = 4 the upper 
Hubbard peak is moved to higher energies. 

Also shown in the graphs are the local particle-hole 
propagators of equation (|36|) (labeled as "free"). These 
show characteristic features of the Hubbard peaks too, 
but the most prominent difference to the fully interact- 
ing susceptibility is the strong suppression of the high 
energy response in the latter. For example, the broad 
excitation continuum in the particle-hole propagator of 
the asymmetric case {U = 4) for energies in the range 
3 < u < 6 is reduced to a very small peak at ^ 3 in the 
fully interacting susceptibility. 

P-^(z^) is just a measure for the phase space volume for 
statistical independent particle-hole excitations, which 
are described by the one-particle DOS pf. The quasipar- 
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tides and -holes at the Fermi level, on the other hand, are 
strongly correlated, leading to an effective suppression of 
the available phase space volume. 

The role of the low energy quasiparticles can be stud- 
ied by comparing the particle-hole propagator and the 
interacting susceptibilities at low energies for the sym- 
metric situation {U = 2, inset). The response via the 
particle-hole propagator for u < Tk shows an increase for 
lower temperatures, which stems from the quasiparticle- 
quasihole excitations within the growing Kondo reso- 
nance. In the fully interacting susceptibility this increase 
is approximately an order of magnitude smaller, clearly 
showing the effect of correlations in the two-particle re- 
sponse. 

Surprisingly, for the asymmetric case {U = 4) this 
trend is reversed and the interacting susceptibility is en- 
hanced over the particle-hole propagator for excitation 
energies smaller than the Kondo scale v < Tk (see in- 
set). This very pronounced low energy response is pro- 
duced by quasiparticle-quasihole excitations in the local 
Fermi liquid phase at low temperatures. The fact that 
these excitations are strongly enhanced in the asymmet- 
ric case compared to the symmetric situation is associ- 
ated with the presence of strong valence fluctuations, as 
was already discussed for the static susceptibility above. 

Even though such a low energy enhancement was 
not reported in NRG calculations for the charge 
susceptibilit y ^^1^^ an inspection of the low energy part of 
the many-body spectrum obtained with the NRGS and 
preliminary NRG-calculations indeed suggest the possi- 
bility of an enhanced charge response for asymmetric sit- 
uations. 

A similar but not as strong temperature dependent 
increase in the dynamic susceptibility for temperatures of 
the order of the Kondo temperature was already found for 
larger orbital degeneracy and U = oo with the NCA 

Therefore we argue, that these findings for the dynamic 
charge susceptibility are in accord with the ones discussed 
in the previous section for the static charge susceptibility. 
The observed increase of the charge response for energies 
smaller than the Kondo scale is indeed physical and due 
to the influence of the valence fluctuations of the low 
temperature Fermi liquid. However, the magnitude of the 
enhancement shown in Figure [13] is arguable, especially 
for the choice of parameters in the present calculation. 

IV. CONCLUSIONS 

We have studied the SIAM within a conserving approx- 
imation, the ENCA, for a variety of model parameters. 
It was shown, that the ENCA constitutes a very accurate 
approximation for the static and dynamic one- and two- 
particle quantities of that model for temperatures down 
to a fraction of the Kondo temperature. It considerably 
improves the Fermi liquid properties and cures shortcom- 
ings of the (S)NCA, like the removal of the divergence of 
limj^^o ^^^"^""ai^) g^-^ 2ero temperature. 



In symmetric situations {2e^ -\-U = 0), the static mag- 
netic and charge susceptibilities were shown to be in ex- 
cellent agreement with the exact Bet he ansatz results. 
This was even true for cases with very small Coulomb 
interaction [/, which could not be expected from the 
beginning, since approximations within direct perturba- 
tion theory with respect to the hybridization usually have 
problems describing the non-interacting case. 

The static magnetic susceptibility is in excellent agree- 
ment with exact Bethe ansatz results in the asymmetric 
situation. This holds also in cases with strong valence 
fluctuations, such as for e-^ = or in the empty- and 
fully occupied orbital regimes. 

However, the static charge susceptibility in the asym- 
metric model is described accurately only in situations, 
where the singly occupied impurity valence state repre- 
sents the unperturbed {Aa = 0) ground state. In addi- 
tion, even though we believe that the qualitative features 
of the charge susceptibility in asymmetric situations are 
captured by the presented calculations, the influence of 
valence fluctuations is probably overestimated for too low 
temperatures. This confirms the expectation, that cross- 
ing diagrams, which are neglected in the ENCA, are es- 
sential for the quantitative description of situations with 
strong valence fluctuations, where the impurity occupa- 
tion is statistically fluctuating. This also is in accord 
with the known pathologies in the one-particle spectral 
function. There, charge and magnetic fluctuations both 
contribute and the overestimation of the charge excita- 
tions at very low temperatures leads to the overshooting 
of the Kondo resonance and the observes spikes in the 
DOS. 

The dynamic magnetic susceptibility is dominated by 
Kondo screening of the impurity spin. The ENCA cor- 
rectly reproduces the temperature and other parameter 
dependencies of the magnetic excitations, and also the 
scaling found in previous studies is obtained. 

The dynamical charge spectrum shows a severe sup- 
pression of high energy excitations due to correlations, 
when compared to the particle-hole propagator, which 
would represent the susceptibility of independent renor- 
malized quasiparticles. Additionally the low energy re- 
sponse for excitation energies smaller than the Kondo 
temperature is also strongly suppressed in the symmet- 
ric case, due to the same correlations between low energy 
quasiparticles. 

In the asymmetric situation the low energy charge re- 
sponse is drastically enhanced and an additional peak 
emerges. This enhancement is attributed to the presence 
of the valence fluctuation fixed point in the asymmetric 
model. Such an enhancement seems quite probable so 
that only the steepness of the increase calculated within 
the ENCA for parameter values chosen is arguable. 

With these findings, the prospects of describing two- 
particle dynamics of lattice systems within the DMFT 
are very promising and results will be presented in a sub- 
sequent publication!^!^ 
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